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Size Evolution of Spheroids in a Hierarchical Universe 
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ABSTRACT 

Unveiling the structural evolution of spheroids, and in particular the origin of the 
tight size-stellar mass relation, has become one of the hottest topics in cosmology in 
the last years and it is still largely debated. To this purpose, we present and discuss 
basic predictions of an updated version of the latest release of the Munich semi- 
analytic hierarchical galaxy formation model that grows bulges via mergers and disc 
instabilities. We find that while spheroids below a characteristic mass M s ~ 10 11 M 
grow their sizes via a mixture of disc instability and mergers, galaxies above it mainly 
evolve via dry mergers. Including gas dissipation in major mergers, efficiently shrinks 
galaxies, especially those with final mass M s < IO^Mq that are the most gas-rich, 
improving the match with different observables. We find that the predicted scatter 
in sizes at fixed stellar mass is still larger than the observed one by up to < 40%. 
Spheroids are, on average, more compact at higher redshifts at fixed stellar mass, 
and at fixed redshift and stellar mass larger galaxies tend to be more starforming. 
More specifically, while for bulge-dominated galaxies the model envisages a nearly 
mass-independent decrease in sizes, the predicted size evolution for intermediate-mass 
galaxies is more complex. The z = 2 progenitors of massive galaxies with M sta r ~ 
(I — 2) x 10 11 Mq and B/T > 0.7 at z = 0, are found to be mostly disc-dominated 
galaxies with a median B/T ~ 0.3, with only ~ 20% remaining bulge-dominated. The 
model also predicts that central spheroids living in more massive haloes tend to have 
larger sizes at fixed stellar mass. Including host halo mass dependence in computing 
velocity dispersions, allows the model to properly reproduce the correlations with 
stellar mass. We also discuss the fundamental plane, the correlations with galaxy age, 
the structural properties of pseudobulges, and the correlations with central black holes. 
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1 INTRODUCTION 

One of the most important and still debated problems in 
Cosmology is the formation and evolution of galaxies. To- 
day we see galaxies having a variety of morphologies, rang- 
ing from less massive, pure stellar discs, to intermediate 
mass bulge plus disc galaxies, to more massive, spheroidal 
systems. The origin of this transition is still for many re- 
spects unclear. More specifically, while angular momentum 
conservation may exp lain many properties of discs (e.g., 
iGovernato et al]|2007l ). the origin of bulges is still largely 
unsolved and debated. Why do some galaxies show bulges 
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while others don't? Or, in other words, what is the origin of 
the gradual conversion from discs to spheroids? 

It is clear that if we want to understand galaxy forma- 
tion we need to observe the high-redshift Universe. Deep 
observations in the last decades or so have however un- 
veiled a full complex zoology of high-redshift (proto) galaxies 
that makes even more puzzling - but also more excit- 
ing - assessing the actual routes chosen by nature to 
build the galaxy populations we observe today. Along 
with starforming discs and dust - enshr ouded galaxies (e.g., 
iBlain et all |2002| ; iMagdis et ail 1201 ll ). deep optical and 
near-infrared surveys have in fact discovered the pres- 
ence of numerous extremely co mpact and passively evolv- 
ing galaxies up to z > 3 (e.g., iTruiillo et all 120061. 120071; 



_ ii) . _ , 

ICimatti et ail boosl : iBuitrago et all 120081 : IChapman et all 
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2008 


; IPranx et al. 2008; 


Saracco et al.l 20081; Tacconi et al. 


2008 


; van der Wei et al. 


2008; van Dokkum et al. 2008; 


Younger et all 12008: 


Bezanson ct al. 2009; Damianov et al. 


2009 ; 1 Williams et"aL 


l20ld; iRvan et al.ll2012l; ISaracco et al. 



20101 ; Ivan Dokkum et alj l2010h . The stellar masses range 



from 10 10 — 10 r2 M@, with a factor of ~ 2 systematic un- 
certainty, and with half-light radiuses within 0.4 — 5 kpc, 
being 2 — 6 times more compact tha n their local coun- 
terparts of similar stella r mass (e.g., Ivan Dokkum et "all 
l20ld ; ISaracco et aT1l201lh . Although several observational 
limit ations may affect h igh-redshift size measurements 
(e.g.. iMancini et"aD l2009h , extremely deep i maging (e.g.. 



Szomoru et al 
2Q111: 



2010h also through lensing fe.g.. lAuger et al.l 



Newton et al.ll201ll ). and the available measurements 



gh velocity dispersions for a subset of these galax- 
ICenarro fc Truiilld 120091; ICappellari et al.1 Bool 



of very hi 

Ivan Dokkum et al.ll2009l ; lvan de Sande et al.ll2011l ). are con- 
firming their extreme compactness. 

Moreover, it is now being established that early-type 
galaxies at higher redshifts are not all compact. At fixed stel- 
lar mass, similar fractions of large and comp act galaxies of 
similar mass co-exist at the sam e epoch (e.g.. lMancini et al.l 
120091 ; IValentinuzzi et alj|2010af), and with a v ariety of bulge- 



to-disc ratios (e.g. 
the spectra o f 62 



van der Wei et all I20TH '). By studying 
early-type galaxies at high redshifts, 



ISaracco et al.l l|201ll ) found that compact galaxies tend to 
have most of their stars formed at z > 5, while larger galax- 
ies at fixed stellar mass are generally younger. Along sim- 
ilar lines and exten ding the anal y sis to other high-redshift 
galaxy populations. iMosleh et all (|201ll ) concluded that the 
structure of galaxies is somewhat correlated to their activity, 
i.e., the sizes of galaxies at a given stellar mass is somewhat 
correlated to its star formation rate level, similarly to what 
is observed in the local Universe. At lower redshifts it has 
been shown that the size-age rel ation at fixed stellar mass is 
similarly shaped for len ticulars (|Shankar fc Bernardill2009l ; 
Ivan der Wei et ah! I2009T ) . i.e., older systems are more com- 
pact, but becomes rat h er flat for bulge-dom i nated galaxies 
JShankar et al l l2010bl ; iBernardi et"afl l2010l : iTruiillo et al.l 
l201ll ). Thus whatever process formed massive ellipticals, it 
must have been fine-tuned to bring all young and old high- 
redshift massive spheroids on the same local size-mass rela- 
tion. 

Understanding the evolutionary link these compact and 
large high-redshift galaxies might have with the variety of 
starforming galaxies at similar redshifts and stellar mass, if 



an open debate (e.g.. Cole et al.l 


2000; [Benson et al. 


2003; 


Granato et al.l |2004|; iMenci et al. 


120041; iBaugh et al.l 


2005; 


Bower et al.1 20061; iGranato et al.ll2006l:|De Lucia et alj2006t 


Menci et al.l 20061; De Lucia & Blaizot 20071; Monaco et al. 


20071; iFan et al. |200SL ISomerville et al.l 120081; iDekel et al. 


2009bl; Khochfar & Silk! 120091; iNeistein & Weinmannl 




Bournaud et al. 201ld; iGonzalez et al.ll201lf). 





According to the standard cosmological paradigm of 
structure formation and evolution, dark matter haloes have 
grown hierarchically, through the continuous merging of 
smaller units into larger systems. In this scenario, galax- 
ies fo rm inside this hi e rarchically growing s ystem of haloes 
(e.g.. ICole et~aT]|200d ; Ibe Lucia et all 120061 ). However, the 
actual role played by mergers (major, minor, wet, and dry), 
in the structural evolution of massive spheroids is still un- 



certain (e.g.. [Hopkins et al.ll201ol and references therein) . 
Some models of g alaxy formation (e.g.. lEggen et al]|l962l ; 
iMerlin et ail 120121 ) envisage that most of the mass in lo- 
cal massive spheroids was formed and assembled in a strong 
and rapid burst of star formation at high redshifts, and the 
remnants evolved almost passively thereafter, without being 
strongly affected by late merging events. 

Galaxy formation models built on top of large N-body 
dark matter numerical simulations or analytic merger trees, 
claim instead that although the stars of the most massive 
spheroids are the oldest being formed at very high redshifts, 
they have assembled a large fraction of their final stellar 
mass only at relatively lat e times via a sequence of minor and 
majo r merger events (e.g., Baugh et al.ll2005l; De Lucia et al.l 
200d; iDe Lucia fc Blaizotl 120071 ; iKhochfar fc Silkl l2006al ; 



Gonzalez et al.l 2011 ). It has long been known that binary 
mergers between discs can indeed produce spheroidal galax- 
ies and also explain many of their structural propert ies (e.g., 
lBarneslll992l ; lHernauistll992l ; lRobertson et alj|2006h . though 
several issues remain to be solved in this basic scenario (e.g., 
iNaab fc Ostrikerl 120091 ). In high-redshift and gas-rich disc 
galaxy mergers, howeve r, gas dissipation inevitably forms 



compact spheroids (e.g., Naab et al]|200d ; iRobertson et al.1 



2006; Hopk ins et al . 2009). Hierarchical models then natu 
rally explain the evolution of the size (and mass) of mas- 
sive compact spheroids as a sequence of "dry" (gas-poor) 
and mainly minor mergers that puff up the outskirts of 
the galaxy leavi ng the central regions of the galax y al- 
most intact (e.g.. INaab et al.ll2007l . 120091 ; ICiottill2009r i. Nu- 
merical simulations have however shed doubts on the co- 
herence with which mergers can bring galaxies along the 
tight structural relations observed in the local Universe (e.g., 
ICiotti fc van Albada|[200ll ; iNipoti et al.ll2009l ) . 

Another class of models explains size evolution of 
early-type galaxies via a quasi-adiabatic expansion phase 
consequent to the blow-out of substant ial amounts of 
mass via quasar and/or stellar feedb ack (|Fan et al.l 120081 ; 
iDamianov et al.ll200d ; IFan et al.ll20ld ) . Initial numerical ex- 
periments to test the latter proposal as a viable expla- 
nation to size evolution have been recently performed by 
iRagone-F igucroa & Granatd |201ll ). 

Bulges could also be formed via in-situ processes that 
are broadly classified as disc instabilities. In unstable self- 
gravitating discs, the instability may drive the formation 
of a bar with mass transferred from the disc into a cen- 
tral bulge (e.g.. lCole et~al]|2000D . The degree of mass trans- 
ferred to the bulge varies from one model to the other. 
Some models consider the instability quite a violent process 
capable of transferring most of the disc in to a bulge and 
also induce a starburst (Bo wer et al.| [2006). Besides bars, 
other types of instabilities could contribute to the forma- 
tion of bulges. Observations from t he deep SINS survey of 
z ~ 2 galaxies (|Genzel et al.l l201ll ; iForster Schreiber et al.1 
l201ll ) have revealed the existence of gravitationally bound 
clumps residing in gas-rich discs. Some hydrodynamic simu- 
lations have confirmed that within the turbulent and gas- 
rich high redshift discs, large gas clumps can indeed be 
formed and migrate via dynamical frictio n to the centre thus 
progressively build a s tellar bulge (e.g.. |Pekel et al.| [2009a: 
iBournaud et al.ll2011al ). 

Given the still debated physical mechanisms behind the 
evolution of spheroids, theoretical work on this hot topic is 
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mandatory. The aim of this work is to explore the full pre- 
dictions of a state-of-the-art semi-analytic model (SAM) of 
galaxy formation that evolves massive spheroids in a hierar- 
chical fashion. More specifically, we will present and discuss 
basic predictions on the size-mass relation and its evolution 
with redshift of relatively massive spheroids. We will mainly 
focus on their size evolution at fixed stellar mass and in dif- 
ferent environments, the structural properties of their pro- 
genitors, but also touch upon several other related issues. 
Our objective is not to prove that the model discussed here 
is the correct one, but rather to lay out the successes and 
failures of a detailed hierarchical model against the wealth 
of data now becoming available from large and deep surveys. 
We will also stress that several outcomes of the model con- 
sidered here are shared by many other hierarchical SAMs, 
making most of our conclusions of particular interest to the 
field of galaxy formation, but also discuss key differences. 

The paper is organized as follows. In Section [2] we in- 
troduce the hierarchical model considered for this work, fo- 
cussing on the features relevant to this work. In the same 
section we present a large sample of well studied early-type 
galaxies taken as term of comparison for model outputs. Sec- 
tion [3] contains a discussion of key issues such as the size 
(and velocity dispersion)-stellar mass relation, its scatter 
and evolution with redshift, and the role of environment. We 
discuss in Section [4] the evolutionary features of spheroids, 
the analogies and differences with other SAMs, and ways 
to further constrain galaxy formation models. We then con- 
clude in Section [5] In the Appendices we will also briefly 
discuss a number of related topics, such as the distinction 
between classical and "pseudobulges" , the connection with 
black holes and with galaxy ages, the fundamental plane and 
its evolution with time. 



2 INITIAL SETTINGS: MODEL AND DATA 
2.1 The Reference Model 

All of the results presented in this work are th e outcome 
of run ning the original numerical source code bv lGuo et al.l 
(|201l] >. i.e., the latest rendition of the semi-analytic model 
(SAM) developed at the Max Planck Institute for Astro- 
physics. As det ailed below, with respect to the original 
IGuo et all (|201ll ) model, we have modified the computation 
of bulge radii exploring a variety of possibilities, and added 
the calculation of the coupled velocity dispersions. Note that 
running the code and producing a new galaxy catalog each 
time, is different from studying the online catalogs as it al- 
lows a self-consistent thorough study of the structural evolu- 
tion of galaxies in the SAM. We stress th at the modifications 
applied to the original IGuo et all l|201ll ) model do not affect 
any other galaxy property except for sizes. Thus the model 
(we checked) maintains the same exact performance with 
respect to the observables ( e.g., the stellar mass function) 
as presented in IGuo et all (^Olll ) . Before discussing bulge 
sizes in detail, we first provide below a brief overview of the 
model. 

The Munich SAM aims at providing a comprehensive 
picture of the evolution of galaxies and their central super- 
massive black holes within the hierarchical structure and 
merging of dark matter haloes and subhaloes within the 



concordance ACDM cosmology. To this purpose, it is imple- 
mented on top of the large, high-resolutio n cosmological N- 
bod y MILLENNIUM I (fSpringel et all 120051 ) and MILLENNIUM 
II (jBovlan-Kolchin et al.ll2009h simulations. Given that in 
this paper we are mainly interested on the structural proper- 
ties of the most massive galaxies in the local Universe, all of 
the results presented here have been obtained by running the 
code on the significantly larger MILLENNIUM I simulation. 
The latter simulation follows the evolution of N — 2160' 3 
dark matter particles of mass 8.6 x I0 8 /["'M©, within a 
comoving box of size 500/i _1 Mpc on a side, from z = 127 
to the present, with cosmological parameters f2 m = 0.25, 
fl b = 0.045, h = 0.73, Ov = 0.75, n = 1, and ag = 0.9 

As comprehensively detailed in IGuo et all (|20 1 if ) , the 
model self-consistently evolves the full population of galaxies 
within the hierarchy of dark matter haloes, adopting a set 
of equations to describe the radiative cooling of gas, the star 
formation, metal enrichment and supernovae feedback, the 
growth and feedback of supermassive black holes, the UV 
background reionization, and the effects of galaxy mergers. 

Particularly relevant to the present paper is the gener- 
ation of galaxy morphology within the model. Collapse of 
baryons within dark matter haloes and conservation of spe- 
cific angular momentum naturally leads to the formation of 
discs. In the Guo et al. model gaseous and stellar discs are 
distinguished, and each component is evolved in time in an 
inside-out fashion, continuously evolving in mass and angu- 
lar momentum due to the progressive addition from cooling 
gas, minor mergers, and gas removal from star formation. It 
is also assumed that both the gas and stellar discs are thin, 
in centrifugal equilibri um, and to have exponential profiles 
(we refer the reader to IGuo et all (|201 lh for full details) . 

One of the primary channels to instead form and evolve 
bulges in the model is via galaxy mergers. In the case of 
a minor merger (Mi/ Mi < 0.3), the disc of the primary 
galaxy survives, and the stars and the gas of the satel- 
lite are added to any pre-existing bulge and to the disc 
of the primary galaxy, respectively. Galaxy major mergers 
(Mi/ Mi ^ 0.3) instead disrupt any stellar disc present and 
produce a spheroidal remnant, which contains all the old 
stars present in the progenitor galaxies and all the new stars 
formed out of the burst triggered by the merger. 

Bulges can be formed even via secular evolution in the 
model. The type of disc instabilities considered in this model 
are secular processes that transfer only the portion of the 
stellar mass necessary to keep the disc marginally stable 
(see details in Section below, and in Guo et al. 2011). This 
way of mod e lling disc instabilities is different from, e.g., 
iBower et al.l (|2006l ) that instead assume the entire mass of 
the disc is transferred to the bulge during the instability, 
with any gas present assumed to undergo a starburst. More 
in general, the present model lacks at the moment any bulge 
formation via strong gas rich disc instabilities and/or clump 
accr etion unde r dyna mical friction. 

IGuo et al.l (|201lh have shown that their model is capa- 
ble of reproducing the size distribution of local discs rea- 
sona bly well, and additional compariso n s can be found in, 
e.g., iFu et all (|2010h : iKauffmann et al.l i|2012l ). In this pa- 
per we will mainly focus on the predicted structural proper- 
ties of massive spheroids and their evolution with redshift, 
and only briefly touch, where relevant, on the evolution of 
discs that will be addressed elsewhere. In particular, our 
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Figure 1. left panel: Predicted median 3D half-mass radius versus stellar mass for different cuts in B/T, as labelled. Right panel: 
median 2D projected half-light effective radius R e as a function of stellar mass for a subsample of SDSS early-type galaxy sample. It 
is clear that the model is at variance with the data, predicting much flatter size-mass relations below a "characteristic mass scale" of 
M a t a r ~ 10 Mq, especially for higher B/T galaxies. This characteristic mass is completely absent in the data. 



primary interest in this work are bulge-dominated galax- 
ies with bulge-to-total ratio B/T > 0.7, a threshold cho- 
sen because, as discussed below, the contribution of pseu- 
dobulges (Appendix [B| in this regime is negligible. This in 
turn allows us to properly discern the actual role played by 
mergers in building their structural properties we observe 
today. Where necessary, we will also devote some attention 
to bulges grown via disc instabilities, though we refer the 
reader to Appendix [Bl fo r some more specific discussion, and 
to the separate work bv lShankar et aJj (2012) for additional 
and complementary analysis of this issue. 



2.2 Computing bulge sizes in the Model 

The bulge 3D half-mass radius Rh of a merger remnant is 
compu ted from energy conservation. Following I Cole et all 
(2000), the model assumes that when two virialized galactic 
systems merge, their Rh is given by 



Efin = £lnt,l + -Bint, 2 + B or b . 



(1) 



The left-hand side of Eq. fTJ is the self-binding energy 
of the remnant, defined as 



B fln 



G(Mi + M 2 f 
Rh 



(2) 



with Rh the half-mass radius of the remnant. 

The terms E- m t,i,2 are the self-binding energies of the 
merging progenitors, and are usually expressed as 



-Bint i 



GMj 

~rT 



(3) 



with Mi the total stellar (including any stars formed during 
the merger) plus cold gas masses. 

The orbital energy B or b is usually expressed in terms of 
the internal energy of the system at the radius of minimal 
separation 



Bo 



/orb GMlM2 

c Ri + R 2 



(4) 



where we initially set / or b = 1 an d c = 0.5 (following 
ICole et al.ll2000l ; lGonzalez et al.ll2009l ). 



Note that the conservation of energy in Eq. |T| is not 
unique and could possibly be expressed in other ways. For 
example, the orbital energy of two galaxies at distance d 
could also be computed in the center of mass system of ref- 
erence as B or b = fiV 2 /2 — GMiM^/d, in terms of the re- 
duced mass /1 = Mi Ma /[Mi + Ma), and relative velocity 
V =|| Vi — V2 II . However, while this expression might have 
the advantage of not being directly dependent on /orb, it 
is still model-dependent. In fact, satellite galaxies stripped 
away of their surrounding subhalo are assigned a surviv- 
ing merging timescale proportional to the actual Chan- 
dras ekhar dynamical fric tion timescale via a fudge factor 
(see lDe Lucia et al.ll201(i |). Thus relative velocities and dis- 
tances would still need to be modelle d according to this 
timescale (see also iNeistein et aTll201ll ). 

Thus, given the inevitable inclusion of some parame- 
ters in the modelling of bulge sizes, we decided for this work 
to stick with Eq. {TJ that requires only one truly free extra 
parameter, / or b /c, and at the same time allows a closer com- 
parison with previous semi-analytic and numerical works in 
the Literature. 

As anticipated in the previous Section, another route to 
form bulges in the model is via secular evolution (Guo et al. 
2011). The adopted criterion for instability is expressed as 
Vniax < \/GMdisc/3-Rdisc, with \4iax the maximum circular 
velocity of the host subhalo, and Mdi sc and Raise the mass 
and exponential scalelength of the disc, respectively. In the 
event of instability, a fraction <5M sta r of the disc stellar mass 
is transferred to a central bulge to restore equilibrium, and 
the size Rb of the newly formed bulge is computed assum- 
ing an exponential profile. If a bulge is already present, the 
size of the bulge is computed via a "merger-type" relation 
between the old and new bulge stellar mass as the one in 
Eq. fTJ, with Mi and Ri the mass and half-mass radius of 
the pre-existing bulge, and M2 and R2 equal to 5M star and 
-Ri,, respectively, and / or b = 2 to take into account that the 
interaction in concentric shells is stronger than in a merger 
(see Guo et al. 2011 for further details). 

After each merger or disc-instability event, we also com- 
pute the velocity dispersion as sociated to each galaxy follow- 
ing the analytical fit given bv lCovington et al.l (|201ll ) 
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2 , GM st 

ar 

a = k- 



Rh 



M H {< R H ) 



M h (< Rh) + /(Afi + Afa) 



(5) 



where we set = 0.15 and / = 0.1. Here Mh{< Rh) is 
the fraction of the subhalo mass associated to the remnant 
within the final half-mass radius. 



2.3 The SDSS sample 

Throughout the paper we will compare model predictions 
with a sample of ~ 25, 000 galaxies from the Sloan Digital 
Sky Survey (SDSS). This is a random subset of the well de- 
fined and complete sample defined in Bernardi et al. (2010). 
This sample ranges between 3 x 10 9 Mq to ~ 10 12 Mq , the 
mass range of interest here. In addition to the photometric 
parameters (e.g., cmodel magnitudes and sizes) presented 
in Bernardi et al. (2010), this subset of ~ 25,000 galaxies 
also p rovides bulge-to-tota l B/T light from a de Vaucoulcurs 
(e.g.. Ide Vaucouleurslll948h plus exponential decomposition 
(see Meert et al. and Vikram et al. in preparation), which 
we will extensively use in this work to select the most bulge 
dominated galaxies to compare with the model. The sam- 
ple is character ized by stellar m asses with a Chabrier Initial 
Mass Function (Chabrier 2003), consistent with the one used 
in the model. 



3 RESULTS 

3.1 A first comparison 

In Figure [1] we show a first comparison between model pre- 
dictions and data from SDSS. The left panel shows the me- 
dian 3D half-mass radius _R.tr versus stellar mass for dif- 
ferent cuts in B/T, as labelled. From here onwards, unless 
otherwise stated, we will compute Rh as the mass-weighted 
average of the half-mass radiuses of the bulge and the disc 
components. We have checked that, especially for the bulge- 
dominated systems of interest here, this is equivalent (at 
the percent level) t o computing a f ull mass profile assuming, 
e.g., an Hernquist ijHernq uist 1990) plus an exponential pro- 
file for the bulge and the disc, respectively. Moreover, simply 
neglecting the disc component in systems with B/T > 0.7, 
the subsample of galaxies this paper focuses on, yields very 
similar results. 

The right panel shows the median effective radius R e 
as a function of stellar mass for the same cuts in B /iQ. The 
model predicts an increasing size with stellar mass, how- 
ever, it is evident that below a "characteristic mass scale" 
of M s tar ~ 10 11 Mq the predicted sizes flatten out, at vari- 
ance with the data that continue to show a steep decline 
down to much lower masses 0- 

It is quite unlikely that such a strong discrepancy can be 



We note that the behaviour of increasing half-mass radius with 
B/T at fixed mass is induced by the fact that the model predicts 
larger bulge and disc sizes with increasing B/T. We verified that 
just the opposite is true for disc- dominated galaxies with B/T< 
0.3 — 0.5, in agreement with observations. 

2 The apparent flattening of the measured size-mass relation at 
very low masses M s tar i; 10 10 Mq is most probably induced by 
contamination of later-type gala xies, and thus not releva nt for the 
present discussion (see details in Bernardi et al. 2011~b|). 



simply explained by some mass/luminosity-dependent con- 
version factor between Rh and R e , as we will also discuss in 
more detail later. We checked that the predicted Rn-Mst^i 
relation shows a very similar flat behaviour at low masses 
already at z ~ 2. The latter implies that the wrong shape of 
the size-mass relation yielded by the model is a consequence 
of some wrong "initial conditions" and not necessarily linked 
with any later galaxy assembly. 



3.2 Including Dissipation 

If the merger is gas-rich, a significant fraction of the initial 
energy -Ediss of the system will be dissipated away, induc- 
ing a more compact remnant. Several groups have studied 
this issue in some detail using high-resolution hy drodynamic 
simu l ations and semi- numerical models (e .g., iNaab et al.l 
l200rj ; ICiotti et al1l2007l ; iHopkins et"al"1l2009l . and references 
therein) . 

Some of these groups have provided basic analytic for- 
mulations that can be included in SAMs to study the impact 
of dissipation on cosmological structure formation of galax- 
ies. One me thod is based on t he co nservation of energy, as 
proposed bv lCovington et alj (|2008l ) 



Efin = -Bint,! + Eint,2 + -E or b + Edis 



(6) 



with the di s sipatio n energy parameterized by 
ICovington et al.1 l|201ll ) in terms of the final energy of 
the remnant as 



-Ediss — 2.75/gas-Efin , 



(7) 



with /g as the ratio between the total mass of cold gas and 
the total cold plus stellar mass (inclusive of the mass formed 
duri ng the burst) of the p rogenitors. 

IHopkins et al.l l|2009h parameterize the decrease in size 
due to dissipation by the simple relation 



Rh [final 



Rh [dissipationless] 

1 + /gas//() 



(8) 



with fo = 0.25, and Rnldissipationless] computed from 
Eq. Q. 

It is clear that Eqs. (JS} and © can significantly reduce 
the sizes predicted in a dissipationless merger. We assume 
that Eqs. (0 and ([8} only hold in major mergers, when 
the gas actually gets into the bulge (see Section |2.2|> . The 
decrease is proportional to the gas fractions in the progen- 
itors. As shown in the right panel of Figure [2l the model 
predicts increasing gas fractions at lower stellar masses and, 
at fixed mass, increasing wi th redshift, in broad agree- 
ment with observations (e. g., Kannappanl 200*41 ; Erb et ah! 
l200fj ; ICatinella etahl l20ld ; iPeeples fc Shankari |20ld and 
references therein). The results shown here refer to galax- 
ies with B/T> 0.5 but they are general to galaxies with 
higher and somewhat lower B/T. Low- mass galaxies can 
easily have most of their baryonic mass still in gaseous 
form at z > 2. However, at any epoch, galaxies above 
M c > 3 x 10 10 Mq tend to have progressively lower gas frac- 
tions down to < 10 — 20%. The mass M c is an interesting 
mass scale several times reported in the Literature to be in- 
dicative of some basic physical process in galaxy evolution 
(feedback from Active Galactic Nuclei?) as several spectral 
and structural properties change when galaxies transition 
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Figure 2. Left panel: The dotted, long-dashed, and solid li nes show the pred i cted m edian size-mass relation for galaxies with B/T > 0.5 
for a model w i th no dissipation, with dissipation following Covington ct al. (2011), and with dissipation following the prescriptions by 
iHopkins et alj (120091 ) , respectively; to vertical dotted line marks the transition above which (dry) mergers are believed to dominate galaxy 
assembly (see text). Right panel: median cold gas fractions predicted in the model as a function of stellar mass at redshifts z = 0, z = 1, 
and z = 2, as labelled. 



above it (e.g., Kauffmann et al.l 2004 


; Shankar et alj|2006t 


iKhochfar & Silk 20091; Bernardi et al. 


l2011alhb. 



This mass- and time-dependent behaviour of / gas can 
then easily explain the decrease in size shown in the left 
panel of Figure The dotted, long-dashed, and solid lines 
show the predicted median size-mass relation for galaxies 
with B/T > 0.5 f or a model with n o dissipation, with dis- 
sipation following ICovington et alj (1201 if), and with diss i- 
pation following the prescriptions bv lHopkins et ail (|2009h , 
respectively. Although a decrease in size is apparent at all 
masses, dissipation is progressively more effective at lower 
masses, proportionally to the increase in gas fractions. In 
both panels for reference the vertical dotted line at M sta r = 
10 11 M© marks the transition above which dissipation does 
not play a significant role in shaping the sizes of galaxies 
because mergers become progressively gas-poorer. It is re- 
markable that dissipation tends to erase the flattening below 
the characteristic mass producing a nearly single power-law 
in agreement with the data. In the following we will use 
Eq. [8] as our reference model with gas dissipation, though 
comparable results are obtained when switching to Eq. [6] 



3.3 A closer comparison to the data 

We now attempt a closer comparison between model predic- 
tions and data measurements by converting 3D half-mass 
radiuses Rh into 2D projected half-light radiuses R e . As- 
suming that light trac es mass we convert Rh to R e using the 
tabulated factors from lPrugniel fc Simienl (|l997l ). The latter 
computed for a system of total mass M, the scaling factors 
S(n), dependent on the Sersic index n (Sersic 1963), con- 
necting the gravitational energy W to their effective radius, 
i.e., \W\ = S{n)GM 2 /R e = GM 2 /R S , with G the gravita- 
tional constant and R g the gravitational radius. Assuming 
the systems are virialized we can approximate R g ~ 2Rh, 
thus having 

Re w 2S(n)R H . (9) 



T ■ 1 I 




NO dissipation 

d d ■ i d i i i i i 

10 10 10" 10 12 
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Figure 3. Effective radius R e as a function of stellar mass 
at 2 = from the SDSS sample of early-type galaxies with 
B/T > 0.7. The contours mark the region of plane containing 
68%, 95%, and 99.7% fraction of the total sample. For complete- 
ness, the solid squares represent the medi an size-mass relation fo r 
the early-type galaxy sample discussed bv lBernardi et al.l ll2011bl ) . 
The dotted, long-dashed, solid, and triple dot-dashed lines repre- 
sent, respectively, the predicted size-mass relations without dissi- 
pation, with dissipation, with dissipation plus / or b = 0, and with 
dissipation plus a fraction of dark matter in the merger (see text 
for details). 

By setting n = 4 in Eq. @ (i.e., S(4) = 0.34 from Table 4 of 
IPrugniel fc Simienl Il997l. we can convert the predicted 3D 
half-mass radiuses into 2D projected half-light radiuse^f). 

Figure [3] shows the z = SDSS effective radius R e 
(assuming de Vaucouleurs plus exponential profiles) as a 

3 We have checked that our conclusions do not significantly de- 
pend on the exact profile chosen for the bulges. For example, 
we find broadly similar results, although somewhat steeper cor- 
relations at the highest stellar masses, when assigning to each 
spheroid a Sersic index acco rding to their lum i nosity following, 
e.g., the empirical relation bv lTerzic fc Gra ham (2005), and then 
converting from Rh to R e using the appropriate S(n). 
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Figure 4. Relative size distribution of galaxies of a given mass and minimum B/T as labelled. The dotted and long-dashed, and solid 
lines are, respectively, for models with no dissipation, with dissipation including and excluding pseudobulges. More specifically, the solid 
lines refer to models with dissipation but include only the subsample of galaxies grown mainly through mergers. Each distribution is 
normalized in a way that the sum of the galaxy fractions in each bin of size equals unity. 



function of stellar mass for galaxies with B/T > 0.7, with 
the contours marking the regions containing 68%, 95%, and 
99.7% of the total sample (the results discussed below do 
not depend on the exact choice of B /T threshold) . For com- 
pleteness, we report in the same Figure with solid squares 
the median size- mass relation fo r the ea rly-type galaxy sam- 
ple discussed bv lBernardi et all l|2011bl ). 

The lines in Figure [3] show the predicted median size- 
mass relations for galaxies selected to have the same mini- 
mum B/T threshold as in the data. The dotted line refers to 
a model with no dissipation. As discussed above, the latter 
model better lines up with the data for galaxies above the 
characteristic mass, but below it the relation inevitably flat- 
tens out to larger sizes. The long-dashed line is the size-mass 
relation for a model that includes dissipation in major merg- 
ers. As expected, the sizes get progressively smaller towards 
lower mass spheroids that have formed out of gas-richer pro- 
genitors. However, spheroids of all masses get shrunk. Thus, 
including dissipation not only steepens the size-mass rela- 
tion, but it also lowers its overall normalization. 

In order to improve the fit to the data with a model 
with dissipation we therefore need to increase the normal- 
ization of the predicted sizes at fixed stellar mass. Fol- 
lowing Eq. @, we have that the size of the remnant is 
Rh cc (Mi + M 2 ) 2 /(Si + E 2 + Sorb + Sdiss), so in order 
to increase the size at each merger event, we need to either 
decrease the denominator, and/or increase the numerator. 
The solid line in Figure [3] is the predicted size-mass rela- 
tion assuming that most of the merger events happen on 
parabolic orbits with null orbital energy, i.e., / or b = 0, a far 
fro m uncommon condition in n umerical simulations. Actu- 
ally, iKhochfar fc BurkertJ (|2006t) studied the orbital param- 
eters of major mergers of cold dark matter halos using a 
high-resolution cosmological simulation finding that almost 
half of all encounters are nearly parabolic. This simple vari- 
ation to the basic model significantly improves the match to 
the data. 

We can however also increase the sizes by assuming that 
the total mass actually taking part in the merger is the sum 



of the baryonic plus a fraction of the dark matter host halo 
mass, i.e., 

Mi = M s tar,i + M coldii + a X Mn(< Ri) (10) 

where Ri is the half-mass radius of the progenitor, and a 
a constant parameterizing the still uncertain effect of adia- 
batic contraction. For each progenitor we take its mass at 
infall and c ompute the frac t ion w ithin the half-mass radius 
assuming a iNavarro et al.l l|l997h profile, an d assigning a 
concen tration following the mean relation bv lBullock et al.1 
l|200ll ). The result is shown with a triple dot-dashed line 
in Figure [3] where we set a = 1. The predicted sizes are 
larger, as expected, though the inclusion of a constant frac- 
tion of dark matter also produces somewhat larger sizes at 
low masses t han actually observed (a similar behaviour was 
discussed bv lGonzalez et ai1l2009T ). 

We conclude that dissipation inevitably shrinks galax- 
ies and thus some additional ingredient must be included 
in the model to reestablish the normalization of the size- 
mass relation. In the following we will use the model with 
/orb = and a = as the reference one, unless otherwise 
stated. We note, however, that including some amount of 
dark matter participating in the merger could still represent 
a viable model if we somehow tune a to properly increase 
with halo/stellar masfl 

3.4 On the scatter around the mean size-mass 
relation 

So far we discussed the median shape of the size-mass re- 
lation. We now turn to the discussion of the dispersion in 
sizes at fixed stellar mass. Figure U shows the relative size 
distribution of galaxies in the mass range 10 < log M sta r < 
10.5 and minimum B/T as labelled, with each distribution 
P(R e |Af s tar) normalized in a way that the sum of the galaxy 
fractions in each size bin equals unity. 

4 We stress here that a model with a ^ in the merger cannot re- 
produce the tilt of the fundamental plane discussed in Section [Cl 
that requires a halo mass-dependent velocity dispersion. 
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Figure 5. Predicted logarithmic 1 — a scatter in sizes at fixed 
stellar at different redshifts, as labelled. Red and cyan lines refer 
to model predictions with and without dissipation, respectively. 
For comparison, also shown the scatter measured in our SDSS 
subsample (thick, dot-dashed line), by Shen et al. (2003, filled 
squares), and by Nair et al. (2011, yellow stripe) for early-type 
galaxies. 



The dotted lines show the P(-R e |M s tar) distributions 
competing to the no-dissipation model. While galaxies with 
B/T> 0.7 (right panel) have a Gaussian-like P(R e \M stav ) 
distribution, galaxies with lower B/T> 0.5 (left panel) tend 
to show a double-Gaussian distribution. The second Gaus- 
sian is characterized by a peak and height a factor of a 
few lower than the first Gaussian, though similar in am- 
plitude. The solid lines show the size distributions predicted 
by the model with dissipation but removing those bulges 
that grew at least 50% of their final size (and a large frac- 
tion of their mass) via disc instability. The latter filtering is 
capable to nearly erase the second lower Gaussian from the 
P(i? e |M st ar) distribution, proving that the latter is not pro- 
duced by mergers but rather by disc instabilities. We find 
that most of the remaining galaxies after the filtering have 
B/T> 0.7. 

At masses higher than > 2 x 10 11 Mq we always find 
single- Gaussian distributions irrespective of the chosen B/T 
threshold. We thus conclude that galaxies of any mass with 
B/T> 0.7 are predicted to possess bulges mainly formed 
through mergers (the long-dashed and solid lines in the right 
panel of Figure [4] coincide). On the other hand, galaxies 
below the characteristics mass M B tend on average to have 
a significant population of pseudobulges with B/T< 0.7. 

A full comparison between model predictions and ob- 
servations is given in Figure \S\ where we plot the predicted 
and observed scatter for galaxies with B/T > 0.7 (so to 
minimize the contribution of pseudobulges in this regime. 
We note, however, that even lower cuts in B/T yield similar 
results, as long as only classical bulges are considered). We 
homogeneously compute the scatter in the data and in the 
model by binning galaxies in stellar mass, build the distri- 
bution in sizes and compute the 68% percentile as repre- 
sentative of its 1 — a dispersion. Predictions for the models 
without and with dissipation are shown with cyan and red 
lines, respectively, at redshifts z = 0, z = 1, and z = 2, as 
labelled. 



The thick dot-dashed line is the scatter derived for our 
SDSS galaxy subsample with the same cut in B/T. We find 
it to be at the constant level of ~ 0.5/ln(10) ~ 0.22 dex 
below M c , steeply decreasing to ~ 0.1 dex, the latter in 
very g ood agreement with the results of Irlvde Sz Bernardil 
(2009) for high-mass early-type galaxies. For completeness, 
the squared point s represent t he sca tter of early- type galax- 
ies calibrated by I Shen et al.l l|2003l ). while the b and marks 
the level of scatter measured, more recently, by iNair et al.l 
(2011). The latter computed the scatter in sizes at fixed stel- 
lar mass for early-type galaxies in rich environments. The 
height of the band represents the uncertainty in scatter they 
claim to find when measuring sizes in different ways. Over- 
all, despite the different techniques and selections adopted 
by the different groups, all measurements agree quite well 
with each other. 

The very first issue to note from Figure [5] is that 
the predicted scatter is comparable to the observed one at 
low masses, becoming progressively larger for galaxies with 
A/star > 3 x 10 10 Mq , the scale M c above which galaxies 
become progressively gas-poorer and their structural evolu- 
tion becomes controlled by dry mergers. The disagreement 
between model predictions and data is mass dependent, and 
contained to be < 40%. This is in line with what claimed 
by several previous semi-analytic, numerical, and observa- 
tional studies that claimed the scatter in size at fixed stellar 
mass to be in disagreement with the observed one, espe- 
cially for bulge dominated galaxies at high stellar masses 
(e.g..lGonzalez et al.ll2009l;lNipoti et al.ll2009l ; IShankar et~aH 
l2010al ; INair et alj|2010l . l201ll ). This discrepancy represents 
a challenge for hierarchical models that needs to be further 
understood. 

When moving to significantly lower cuts in B/T the 
significant contamination from pseudobulges tends to fur- 
ther increase the disagreement with observations. We further 
note that the level of predicted scatter is only marginally 
dependent on dissipation. In fact, galaxies with stellar mass 
Af s t a r > 3 x 10 10 Mq have comparable levels of scatter in 
the two models, as expected given that the role of (dry) 
mergers in the size evolution becomes progressively more 
important at higher masses (see Section 14. ip . The latter 
feature is in broad agreement with some previous studies, 
though the role of dissipation in det ermining the final scatter 
of sizes was more emphasized (e.g.. iKhochfar fc Silklkooebl ; 
ICovington et al]|201ll ). Also, the scatter does not strongly 
depend on redshift, especially for more massive galaxies, 
again those with M star >3x 10 10 Mq. 



3.5 Size-evolution with redshift 

As discussed in Section [T] massive and passive spheroids at 
high redshifts appear more compact with respect to their 
local counterparts in SDSS. This section is dedicated to un- 
derstand the degree of size evolution at fixed stellar mass 
predicted by the model and discuss it in light of the avail- 
able data. 

Figure [7] shows the predicted size-mass relation from 
our reference model, in terms of the projected radius R e , at 
different redshifts, as labelled, and compared to SDSS local 
values. Here at all redshifts we select only bulge-dominated 
galaxies with B/T > 0.7. We find a progressive decrease 
of sizes at all masses. Galaxies with stellar masses above 



© 2012 RAS, MNRAS 000,[L>[2T] 



Size Evolution 9 



10.0 r 



u 

Q. 



or 




Figure 6. Predicted median size-mass relation for bulge- 
dominated galaxies with B/T> 0.7, in terms of the projected 
radius R e , at different redshifts, as labelled, and compared to 
SDSS local data (coloured regions). There is a progressive de- 
crease in sizes of up to a factor of ~ 2 at z = 3, similar at all 
stellar masses. 



M star > M c = 3 x 10 10 M seem to experience a similar de- 
gree of evolution, i.e., a progressive average decrease in size 
U P to < 3 at z < 3. This type of nearly mass-independent 
evolution implies that the slope of the size-mass relation is 
predicted to be almost constant at all redshifts, at least for 
bulge-dominated systems with B/T > 0.7 with M ata r > M c . 

Figure [7] presents a more specific study of the redshift 
evolution in the median half-mass radius Rh, normalized to 
the median local value, for two different intervals of bulge- 
to-total ratios, i.e., B/T< 0.4 and B/T> 0.7, limits chosen 
to select statistically significant disc- and bulge-dominated 
galaxy samples, respectively. Each point in the Figure repre- 
sents normalized median values with their associated error 
bars. The upper panel shows the evolution competing to 
lower mass galaxies with 10 10 < A/ Star /M < 3 x 10 10 , the 
middle panel for galaxies with 3 x 10 10 < M 8ta r/M Q < 10 11 , 
while the lower panel refers to higher mass galaxies with 
10 11 < M star /M < 3 x 10 11 . Globally, for all galaxies, even 
for the bulge-dominated ones (solid lines), we do not find 
strong dependence of size evolution on gas dissipation (red 
and cyan lines refer to model outputs with and without dis- 
sipation, respectively). We have also checked that the degree 
of evolution in bulge sizes does not significantly depend on 
the amount of orbital energy included in the model (Eq[T]). 
This is not unexpected given that the degree of bulge size 
evolution in hierarchical models is mainly governed by the 
number and type of mergers (Section 14. ip . In other words, 
dissipation mainly acts in deciding how compact spheroids 
appear after the initial gas-rich major merger event, leaving 
the degree of evolution, controlled by other processes, not 
significantly perturbed. 

From Figure [7] it is evident that the model, as antici- 
pated in Section 12.11 can manage to reproduce the moder- 
ate redshift evolution in the half-mass stellar radius of disk- 
dominated galaxies (long-dashed lines), which are observed 
to d ecrease in size by a fa ctor of ~ 2 up to z — 2 (see, 
e.g.. ISomerville et al.ll2008l . and references therein). On av- 
erage, however, we do not find the empirical trend for which 





1 


2 : 




1 


o : 














ii 

IM 


n 
\j 


H - 

(J 


I 






ct: 
\ 





6 : 








N 






x 





4 7 


rv 











2 : 







o ■ 




1 


2 : 




1 


o : 








o~ 






II 

N 


n 


H - 

u 


I 






rr 

\ 





6 : 








N 






I 





4 ; 


rr 











2 : 







o : 




1 


2 ; 




1 


o : 








o" 






II 

M 


n 

\J 


8 - 


X 






OT 
\ 





6 : 








N 






X 





4 : 


rr 











2 : 







o ■ 



10.0<Loq M STAR <10.5 0.0<B/T<0.4 - -- 
0.7<B/T<1.0 



dissipation 
NO dissipation 




Ryan et al. 



■ Newman et al. 

' ' ' ' I i i i i I ' ' ' ' I 
10.5<Log M STJR <1 1.0 




-+- 



1 1.0<Log M STAR <1 1.5 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 
z 

Figure 7. Predicted median redshift evolution (normalized to the 
median value at z = 0) for the total half-mass radius of galaxies of 
different stellar mass as labelled and B/T as labelled. The upper, 
middle and lower panels refer to, respectively, the size evolution 
of low, intermediate, and high stellar mass galaxies divided into 
two intervals of B/T, as labelled. For bulge-dominated galaxies, 
the redshift evolution does not depend much on stellar mass or the 
degree of gas dissipation. Massive galaxies tend to evolve slower 
than what the data suggest (black, solid and long-dashed lines). 



bulge-dominated galax ies tend to grow faste r than disc- 
dominated ones (e.g.. ISomerville et all I2008I : iRvan et al] 
|2012| ; iHuertas-Companv et al.ll2012T ). If anything, we start 
seeing a flatter size evolution in disc-dominated galaxies only 
for galaxies above > 10 Mq and at z > 1. 

More generally, the model does not predict a very 
strong size evolution for bulge-dominated galaxies at fixed 
stellar mass, which might not be in full agreement with 
the data. From the observational point of view in fact, 
despite the complexities in th e diffe r ent selection pro- 
(e.g.. | van Dokkum et all l2010l: Ivan der Wei et al] 



l201ll ; IHuertas-Companv et al]|2012i ; iNewman et al.ll2012l ). 

passive early-type galaxies tend to show a stronger size evo- 
lution at high stellar masses, decreasing by up to a factor of 
< 5 (e.g, Ivan Dokkum et all 120101 ; IHuertas-Companv et al] 



© 2012 RAS, MNRAS 000. [11-1211 



10 F. Shankar et al. 




10" 10 12 0.0 0.5 1.0 1.5 2.0 

M STAR [M e ] z 



Figure 8. Left panel: Predicted median size-mass relation for all central {dot-dashed line) and Type 2 satellite galaxies {long-dashed line; 
with no restriction in host halo mass), compared to SDSS data. A systematic difference of ~ 30% with centrals being larger than satellites, 
is apparent at all masses. The dotted and solid lines refer to the size-mass relation of centrals in haloes with mass 10 13 < Mh/Mq < 10 14 
(Groups) and Mh > 10 14 Mq (BCGs), respectively. Central galaxies residing in more massive haloes tend be larger. Right panel: median 
fractional size evolution for BCGs, centrals in groups, and Type 2 satellites. BCGs tend to have a much faster evolution than all other 
galaxies of similar mass. 



|2012| ). We report in each panel of Figure [7] the analytic fit of 
normalized size evoluti o n at fixed stellar mass empirically 
derived by iRvan et al.l |2012, black solid lines), based on 
direct deep H imaging and data from the Literature. The 
model predictions are in good agreement with the empirical 
fit for masses below ~ 3 x 10 10 Mq , where the inferred size 
evolution is weaker, but tend to progressively depart from 
the data at higher masses (we stress that this comparison is 
still at a qualitative level as the fits to observations have not 
been derived for a homogeneous sample of galaxies selected 
to have the same B/T as in the model). 



Similar results on the possible inefficiency of merg- 
ers in puffing up massive galaxies ha ve now been 
claimed also by several independent works Cimatti et al.l 
ll2012t);lHuertas-Companv et all (|2012l ); lNipoti et alj (|2012f ). 
iHuertas-Companv et al.f l 2012? ) more recently claimed evi- 
dence for a weaker size evolution at intermediate masses 
3 x 10 10 < M st ar/M Q < 10 11 , in better agreement with 
model predictions, but tend to confirm the strong drop 
in sizes for the highest stellar mass bins. For complete- 
ness, we also report in the bottom p a nel of Figure [7J the 
fit recently inferred by iNewman et all (|2012t ) (long-dashed 
line) that better lines up with the model predictions, at 
least at z < 1. We caution, however, that due to their 
broader selections, their sample may not be restricted to 
only early-type, bulge-domina ted galaxies (see discussion in 
IHuertas-Companv et al.ll2012l . and references therein). 



In our study of size evolution at fixed stellar mass we 
also tried to separate spheroidal galaxies that have mainly 
grown their bulges via mergers from those that mainly grew 
their bulges via disc instabilities. We found tentative evi- 
dence for pseudobulges to evolve slower in sizes with respect 
to classical bulges of similar stellar mass, but the statistics in 
some bins is poor and the systematic difference is confined at 
the < 20% level. SAMs that adopt stronger disc instabilities 
could provide different conclusions in this respect. 



3.6 Role of Environment 

Not all early-type galaxies may follow the same size-mass re- 
lation. Environment, or simply the special location of galax- 
ies at their formation epoch, might induce different evolu- 
tion at later times in galaxies of similar mass. For exam- 
ple, if mergers dominates the structural growth of galax- 
ies, at least at lower redshifts, then galaxies in denser en- 
vironments where mergers are more efficient might appear 
larger at fixed stellar mass. From the observational point 
of view this is still debated. At high redshifts, while some 
groups find clear evidence for larger ga laxies in denser en- 
vironments, at fixed stellar mass (e.g., ICooper et al. I l2012l ; 
IPapovich et a,l.|l2012l). other don't or c l aim some stellar mass 
depe ndence ( Huertas-Companv et alj |2012| ; iRaichoor et al.l 
l2012f ). In the local Universe, when selecting galaxies of a 
given stellar mass in the field and in overdense regions such 
as Clusters, two main issues have emerged recently. Galaxies 
in Clusters less massive than ~4x 10 11 Mq tend to appear 
smaller a t fixed stellar mass than t heir local counterparts in 
the field (IValentinuzzi et aill201 0a) , and non-central cluster 
galaxies might have had a slower or even negligible evolu- 
tion d own to z = (IValentinuzzi et al1l2010bl ; ISaracco et alj 
2010). The observational evidence though is still sparse or 
not secure. IWeinmann et all (|2009l ) investigated size distri- 
butions, among other properties in the SDSS Data Release 
4 group catalogue of I Yang et alj (|2007l ). finding no clear 
difference in the sizes of early-type centrals and satellites. 
Even the degree of evolution for the b rightest cluster g alaxies 
(BCGs), is not yet well understood. iBernardil (|2009r ) found 
that BCGs have evolved by 50% in size in the last few Gyrs, 
while IStott et al. I (l201ll) c laim a milder evolution of > 30% 
since z = l. lAscaso et al.1 l|201ll ) also claim evidence for sig- 
nificant size evolution though not in the light profile (the 
measured Sersic index is nearly invariant with time). 

Given the non-trivial impact that environment might 
induce on galaxy size evolution, it is thus mandatory to 
study what the model predictions are with respect to this 
important issue. We recall that the lGuo et al l (|201ll ) model 
follows in great detail the fate of the gas and stellar compo- 
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nents of galaxies becoming satellites in larger dark matter 
haloes. Satellites suffer tidal and ram pressure stripping that 
can remove a large part of their gas reservoir. In particular, 
the model further assumes that when the host subhalo is 
completely disrupted during its journey within the largest 
halo, the ga laxy's stellar com ponent starts also being dis- 
rupted (see iGuo et all l201ll for details). The stars which 
are stripped away can then later become part of the central 
galaxy of the parent dark matter halo hosting the satellite 
galaxy. Given all of these physical prescriptions, it is nat- 
ural to expect some structural differences between central 
and satellite galaxies in the code. 

The left panel of Figure [8] shows the predicted size-mass 
relation for central and satellite galaxies (dot-dashed and 
long-dashed lines, respectively). Here satellites are only the 
ones defined to be "Type 2" in the code, i.e., the ones that 
have completely lost their associated subhalo due to disrup- 
tion. We find a relatively small, although systematic, differ- 
ence with satellites being smaller by < 30% with respect to 
centrals of the same stellar mass, s omew hat in between the 
finding s of IWeinmann et all (2009) and IValentinuzzi et all 
|2010al ). The dotted and solid lines refer to the size-mass 
relation of centrals in haloes with mass 10 13 < M H /M© < 
10 14 and Mh > 1O 14 M0, respectively. Central galaxies re- 
siding in more massive haloes tend to be larger mainly be- 
cause they have undergone a larger number of mergers over 
cosmic time. 

We have also analyzed size evolution for galaxies liv- 
ing in different environments. The right panel of Figure [8] 
shows the median fractional size evolution for galaxies at 
fixed stellar mass and environment. We consider galaxies 
having stellar masses above > 10 11 M and B/T > 0.7. 
We consider BCGs, identified as centrals in haloes with 
Mh > 10 14 Mq (solid line), centrals in galaxy groups with 
10 13 < Mh/Mq < 10 14 (dotted line), and Type 2 satellites 
galaxies with no restriction in host virial mass (long-dashed 
line). The model predicts that, at fixed stellar mass, galaxies 
residing in progressively more massive haloes have a propor- 
tionally stronger size evolution evolution, mainly induced by 
the larger number of mergers. In particular, BCGs are pre- 
dicted to increase in size by > 50% at z < 1, a degree of evo - 
luti on in between t he one calibrated by IStott et al.l (|201lh 
and lBernardl (|2009h . 



3.7 Additional constraints from velocity 
dispersions 

Several additional clues on the hierarchical evolution of 
spheroidal galaxies can be obtained when considering ve- 
locity dispersions and the virial mass of galaxies. In this 
section we will discuss the physical implications that can be 
derived from the comparison of model predictions with a va- 
riety of key observables that include velocity dispersion. To 
this purpose, unless otherwise stated, we will only consider 
here the subsample of pure elliptical galaxies, i.e., those with 
B/T 0.9 for which a clear definition of velocity dispersion 
is possible and less biased, both observationally and theo- 
retically, given the near absence of a stellar disc (although 
we note that a less extreme cut in B/T does not alter the 
conclusions below). 

All 3D velocity dispersions are computed following 
Eq. ((51 and then converted, consistently with what we dis- 
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Figure 10. Predicted median a - V v i T relation at different red- 
shifts, as labelled, for a model with a dark matter mass-dependent 
a and galaxies with B/T > 0.9. The grey strip e indicates the ve - 
locity dispersion-circular velocity correlation bv lBaes et "all ||2003| 1 
for early-type galaxies with circular velocity converted to veloc- 
ity at the virial radi us using the velocity-dependent correction of 
iDutton et all l|2010l l. 



cussed f or sizes (Section 13. 3p, to l ine-of-sight ID ct(1D) us- 
ing the IPrueniel fc Simienl fl997) Sain) coefficients, i.e., 
a(lD) = {3S K {n)} 1/2 a(3D), where we set n = 4 (the ef- 
fect of the latter correction is however relatively small and 
does not minimally alter our conclusions). 

The very first correlation with velocity dispersion that 
is usually studied is the one between luminosity /stellar 
mass and velocity disp ersion, the Faber- Jackson relation 
|Faber fc Jacksonlll976l ). Our results are presented in Fig- 
ure [9] where the predicted M st ar-o" relation is plotted at the 
redshifts z — 0, 0.5, 1, 2, 3 with a solid, dashed, dotted, long- 
dashed, and dot-dashed line, respectively, as labelled. The 
grey band is the measured M s tar-cr relation from our sub- 
sample of SDSS galaxies with B/T > 0.9. The left panel of 
Figure [9] shows predictions for a model with a computed via 
Eq. ([5]), that includes dependence on the host dark matter 
halo, while in the right panel a simply scales inversely with 
half-mass radius, oc M sta , T /RH- It is clear that the model 
with a computed via Eq. ((5)l provides a much better match 
to the data, with stellar mass properly increasing with in- 
creasing a. Neglecting any mass dependence in a inevitably 
produces a flattening at high masses where the model pre- 
dicts a quasi-linear dependence between half-mass radius 
and stellar mass (cfr., e.g., Figure [3}- It is also interest- 
ing to note that neglecting any halo mass dependence in a 
produces a much stronger evolution in velocity dispersion at 
fixed stellar mass, while a model with halo mass dependence 
contains the evolution in a to < 30% in good agreement with 
direc t observations (e.g..lBernardj|2009;ICenarro fc Truiillol 



I2OO9I ; ICappellari et al.ll2009l ; Ivan de Sande et al.ll201ll ). 

Several groups have proven that a clear correlation 
exists between velocity dispers ion and circular velocity at 
the outer optical radius ( e.g.. |Ferraresel [2002; Bacs et all 



|2003| ; iPizzella et ail 120051 ; IChael 12011^ This is expected 
from basic dark matter theory (e.g., lLoeb fc Peebles! [20031: 



2006; 



Cirasuolo et al. l l2005l : IShankar et al 
2009) as during the early fast-collapse phase of a 



Lapi fc Cavalierel 
halo, 
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Figure 9. Left: Predicted M star -cr relation at different redshifts, as labelled, for a model a dark matter-dependent <r as suggested by 
numerical simulations {left panel) and a model with no dependence on dark matter {right panel). Only galaxies with B/T > 0.9 are 
considered. The predictions in both panels are compared to the subset of SDSS galaxies with same B/T cut {grey areas). 



its potential well is established and a dynamical link be- 
tween baryon velocity at the center and halo circular ve- 
locity should be established. Several high- resolution simula- 
tions have conf i rmed this behaviour (e.g ., IZhao et al]|2003l ; 
iDiemand et"afl 120071 ; IWang et all 120111) . We plots in Fig- 
ure[lO]the predicted median correlation between ID velocity 
dispersion and v irial velocity V v i r of the halo at in fall (see 
iGuo et alJlioill ), at different redshifts, as labelled, for the 
model with a dark matter mass-dependent a and galaxies 
with B/T > 0.9. The grey stripes in Figure [10] indicate the 
correlation between vel o city d ispersion and circular velocity 
inferred by iBaes et al l l|2003l ) for early-type galaxies with 
circular velocity converted to velocity at the virial radius 
us ing the empir i cally-d erived velocity-dependent correction 
of iDutton et alJ (|2010l ) for early-type galaxies (their Eq. 3). 
A good agreement is found in the local Universe. The model 
then predicts some evolution in the zero point of this rela- 
tion with galaxies at fixed velocity dispersion being mapped 
into haloes with higher virial v elocity but the correlation is 
preserved despite mergers (e.g., Bovlan-Kolchin et al .1 [20051 : 
iRobertson et alJl200d : ICiotti et al.ll2007h . We also find (not 
shown in the Figure) that the scatter in velocity dispersion 
at fixed Vvit increa ses with increasing redshift, in line with 
some observations (|Courteau et al ] |2007l ; lHo|[2007l ). 



4 DISCUSSION 

4.1 Evolutionary features 

Having described here and in the previous sections global 
structural properties of spheroids of different masses and 
B/T, we now attempt to sketch a more comprehensive pic- 
ture of spheroid evolution. To this purpose, we select sub- 
samples of 100 galaxies of a given stellar mass and B/T at 
2 = and trace back in time the most massive progenitor 
of each galaxy and record its properties. 

The result is given in Figure [11] which shows the ex- 
pected median redshift evolution of several properties char- 
acterizing spheroid progenitors. We show results for three 
subset of galaxies classified at 2 — by having 10 < 
logMstar < 10.3 (long-dashed lines), 11 < logM star < 11.3 



(solid lines), and logMstar > 11.5 (dotted lines). All galax- 
ies have B/T > 0.7 at 2 = 0. The model predicts that 
the progenitors' half-mass radiuses shrink when moving to 
higher redshifts (top panel; note that we are here plotting 
progenitors identified from the merger trees, while before 
we always considered different galaxies of the same stellar 
mass at different epochs). Noticeably, all galaxies that end 
up being large spheroids in the local Universe are predicted 
to share, on average, quite similar size evolutions, at least 
at 2 < 1.5 — 2, though different morphologies at higher red- 
shifts. 

We find that most of the progenitors of bulge-dominated 
galaxies at high redshifts are disc-dominated (middle panel). 
More specifically, galaxies that today lie below the character- 
istic mass of 10 Mq , are found to turn into disc-dominated 
systems at 2 ~ 1.5 — 2, with B/T < 0.2, and extremely 
gas-rich (Figure [2]). Even more massive galaxies with stellar 
mass at z = in the range 10" < M st ar/M < 2 X 10 11 Mq 
rapidly turn into discs with a median B/T ~ 0.3. In fact, 
we found that only 20% of the galaxies in this mass regime 
remain bulge-dominated with B/T > 0.7 at higher redshifts, 
with the majority turning into disc systems with B/T ~ 0.5 
by z > 1.5 - 2. 

The model also predicts that at all redshifts the most 
massive galaxies are the least starforming galaxies, up to an 
order of magnitude at z < 2 (bottom panel). The model 
therefore suggests that the observe d high redshift gas-rich, 
starforming, and clumpy discs (e.g.. iForster Schreiber et al.l 
Eon]), may be good candidates for being progenitors 
of today's intermediate massive, early-type galaxies with 
Af s tar (2 = 0) < 10 11 Mq, while the compact and red galax- 
ies should end up being the most massive ellipticals we ob- 
serve in the local Universe. Full exploration of high and low 
number densities of galaxies of a given property will help to 
further constrain the model (see Section f4.3p . 

We also note that beyond z > 2 galaxi es are observed 
to sh ow a flattening of their SSFR (e.g., Gonzalez et al. 
1201 ll ). while the predicted SSFR still continues to raise be- 
yond 2 = 2 (bottom panel). Solving the discrepancy between 
model and data is beyond the scope of the present paper. 
We note however that it is a common feature of many galaxy 
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Figure 12. Fractional average mass {left) and size (right) growth of the subsets of galaxies discussed in Figure ITT1 Long-dashed, solid, 
and dotted lines refer, respectively, to the fractional cumulative growth down to z = as a function of final stellar mass, experienced by 
galaxies via major mergers, minor mergers, or any other mechanism (disc instability, star formation, etc.). 



formation m odels, and it has been re centl y pointed out and 
discus sed by IKhochfar fcljiikl \ 201 ll ) and IWeinmann et all 
|201ll ). 



Figure [121 shows the fractional average mass (left) and 
size (right) growth of the subsets of galaxies discussed in 
Figure [ll] Long-dashed, solid, and dotted lines refer, respec- 
tively, to the fractional cumulative growth down to z — 
as a function of final stellar mass, experienced by galax- 
ies via major mergers, minor mergers, or any other mecha- 
nism (such as disc instability and/or in loco star formation). 
We find that the massive spheroids considered here mainly 
grow through mergers in the model. Above M sta r > 10 11 Mq 
these galaxies grow more than > 50% of their final stel- 
lar mass and size via minor mergers, while major mergers 
dominate the growth at lower masses and become progres- 
sively less important at higher masses fe.g.. lKhochfar fc Silkl 
2009). This is expected given that the median accretion his- 
tory of their typical host dark matter haloes, in the range 
10 12 - 10 13 h^M®, are found in high-resolution numerical 
simulations to be dominated by merge rs with satellites tha t 
are ~ 10% of the final halo mass (e.g.. IStewart et al.| [2008). 



We also find that most of the minor mergers, especially 
in the most massive galaxies, are dry, i.e., have a (cold) gas 
mass fraction in the progenitors that is lower than 0.15. Mi- 
nor dry mergers can roug hly preserve the pr ojections of the 
fundamental plane (e.g., ICiotti et alj 120091 . and references 
there i n), and galactic cent ral densities (e.g., ICimatti et al.l 
2008; Bez anson et al.ll2009l . and references therein), though 
they may also increase the scatter in the scaling relations 
(e.g., iNipoti et alj |200SL and references therein). The find- 
ings of Figure [12] is consistent with an inside-out evolution- 
ary scenario, where stellar matter is continuously added to 
the outskirts of the compact h i gh-redshift galax ies as time 
goes on (e.g. jNaab et alj|2009l ; lOser et alJl201Ch . More de- 
tailed comparisons with metallicity, age and colour gradients 
are needed to set this model on firmer foots (see discussion 
in Section l4~3ll . 



4.2 Comparison with other models 

While all galaxy formation models have discussed predic- 
tions for stellar mass distributions, just a handful have taken 
a step further to also consider structural properties, espe- 
cially for early- type galaxies. 

IKhochfar fc Silkl §006b) within the context of a full 
semi-analytic model, emphasized the role of gas dissipation 
in mergers in forming compact massive spheroids at high 
redshifts, deriving progenitors' gas fractions in good agree- 
ment with the ones discussed here. 

lAlmeida et al.1 (|2007l . 120081 ), following ICole et all 
(2000), studied the galaxy size- mass and global fundamental 
plane relations predicted by two significantly di fferent rendi- 
tions o f the GALFORM semi-an alytic model bv lBaugh et al.l 
l|2005l) and iBower et al.l l|2006l ). After varying most of the 
parameters relevant for determining bulge sizes, they con- 
cluded that both models fail to reproduce the sizes of bright 
early-type galaxies, though they noted that a better match 
to the data was ach i eved with no adiabatic contraction (e.g., 
iTissera et all |2010| ; ICovington et all 1201 ll). as assumed in 
this work . On a similar note, Gonzalez et al.l l|2009l ). fol- 
lowing on lAlmeida et ah! (120071 ). further studied the galaxy 
size-mass relations in the GALFORM models varying other 
parameters, such as orbital energy in the merger, but still 
finding significant disagreement with the data. 

IShankar etaD l|2010bl ) and IShankar etahl (|2010al ) 
showed that the size-mass relation at z = predicted by 
the Bower et al. model is much flatter than the observed 
one due to too large low-mass galaxies with stellar mass 
Mjtar < 10 11 Mq, similarly to what found here. They sug- 
gested that the latter problem may be linked to the ini- 
tial conditions, given that large and low-mass galaxies are 
present at all epochs in the model, in line with what dis- 
cussed here. 

We have also adapted the recipes for bulge size growth 
via mergers and disc instability discusse d in Secti on l2.2l to a 
previo us version of the Munich co de by Dc Luci a fc BlaizotJ 
(2007). The latter differs from the lGuo et"aH l|201ll ) version 
in several respects, from physical recipes to values of the 
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Figure 11. Evolution with redshift of the average properties 
of a subset of 100 galaxies derived by following the most mas- 
sive progenitor back in time along its merger tree. We show re- 
sults for three subsets of galaxies classified by having at z = 
10 < logA/ star < 10.3 (long-dashed lines), 11 < logM sta r < H-3 
(solid lines), and logAf s t ar > 11.5 (dotted lines). All galaxies have 
B/T > 0.7 at z = 0. The top panel shows the median fractional 
size evolution, the middle panel the average bulge-to-total ratios, 
the bottom panel the specific star formation rate. 



best-fit parameters. We have checked, however, that most of 
the results discusse d here remain globally similar. 

iMonaco et al] l)2007l ) implemented the recipes from 
ICole et all (|2000h for size growth in their MORGANA model 
(with no dissipation) finding acceptable agreement with the 
local size-mass relation for early-type galaxies, though with 
a scatter larger than the observed one (see their Figure 16). 
Understanding the success of the latter model with respect 
to the previous mentioned ones relies on pinning down the 
differences in the physical inputs of the MORGANA model 
and the impact they have on size evolution. 

Other works aimed at studying the global structural 
properties of local early-typ e galaxies throug h semi -analytic 
techniques was purs ued by ICirasuolo et al] (|2005l ). within 
the framework of the lGranato et al.l (|2004r ) model for the co- 
evolution of super-massive black holes and their host mas- 
sive spheroids. Cirasuolo et al. showed that by tightening 
velocity dispersion of spheroids to t he virial velocity at th e 
epoch of their formation (see also lLoeb fc Pee bles 2003), 
the local early-type velocity dispersion function (ISheth et al.l 
2003; IShankar et alJliooj ; iBernardi et~aHl2010() and Gaus- 
sian dispersion in sizes was fully recovered. Their study is 
quite intriguing as the match to the photometric and dy- 
namical properties of local ellipticals only relies on galaxy 



properties at virialization epoch and minimizes the role of 
later merger events. 

A more refined theory of structural ev olution of 
spher o ids be sides mergers, has been presented by (|Fan et al.l 
I2008L 120101) . As anticipated in Section [T] this class 
of models explains size evolution of spheroids via 
expansion consequent to the blow-out of substantial 
amounts of mass via quasar and/or stell ar feedback 
llFan et al. 1 120081 . |2010| ; iDamianov et al. I l2009l ). Numerical 
support to the latter models wa s recently provided by 
iRagone-Figueroa fc Granatol l|201ll) who showed that even 
in the presence of dark matter, baryons can indeed expand 
by a factor of a few consequent to substantial mass losses. 
They also pointed out that the puffing up via expansion may 
be too rapid with respect to the old ages measured for the 
compact high-z early-type galaxies. Understanding the ac- 
tual role played by expansion versus mergers is beyond the 
scope of the present work. However, we discussed in Sec- 
tion [33] that the strong size evolution for the most massive 
ellipt i cals claimed by some observational groups |Rvan et all 
120121 ; iHuertas-Companv et ail |2012| ; iNewman et al.l |2012| ) 
can hardly be reproduced in our only-merger model, and 
possibly some extra-expansion (a factor of ~ 2) in the first 
phases of evolution might help. 

Relevant numerical work has been pursued in the last 
years to explore the size evolution of spheroids i n a full cos- 
molog ical c ontext. We r ecall here the work by iNaab et al.l 
(|2009h and lOser et al.l (|20ld . see also Scannapieco et al. 
2011), who developed high resolution hydrodynamical cos- 
mological simulations of massive spheroidal galaxies. They 
particularly emphasized that galaxies above the character- 
istic mass of A/ S tar > 10 Mq can accrete via minor merger 
about 80% of their final stellar mass, in agreement with what 
found here (Figure [12]). 



4.3 General issues and additional constraints 

In this work we emphasized the role of gas dissipation as a 
viable mechanism to shrink bulge sizes, especially in lower 
mass galaxies, thus improving the match to a variety of dif- 
ferent observables. Gas dissipation is indeed a natural out- 
come of gas-rich mergers and thus should to be properly 
included in complete models of galaxy formation. Neverthe- 
less, we cannot rule out that part of the discrepancy between 
model predictions and data, especially regarding the match 
with the size mass relation (Section 13. 2p . could also be in- 
duced by some other wrong model inputs. For example, it 
has been recognized that this SAM, like several others, over- 
produces the stellar mass function at high redshifts and low 
stellar masses fe.g.. Irlenriques et al.|[2012l ). This in turn im- 
plies more massive and larger galaxies in lower mass haloes 
with respect to what expect ed from, e.g., cumu lative abun- 
dance matching arguments (jMoster et al.ll2012l ). thus possi- 
bly contributing to the flattening at low masses in the size 
mass relation (Section [3TTJ . Besides the actual performance 
of the model in properly predicting the size-mass relation, 
we nevertheless stress that gas dissipation can improve the 
match to several other observables, as discussed above and 
further in the Appendices (e.g., size-age relation, correla- 
tions with velocity dispersion, etc.). 

In this work we have focused our attention on the most 
relevant median scaling relations among structural proper- 
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ties of bulged galaxies. It is clear, however, that several other 
properties are equally important to pin down and better con- 
strain the viable models. 

First of all, it is fundamental for a model not only to 
produce the correct structure but also the correct number 
of galaxies of a given type. Galaxy formation models have 
most seriously investigated the match to the stellar mass 
function, finding good agreement adopting different physi- 
cal prescriptions (e. g., Cole et al]l2000l; iBenson et al.ll2003l ; 



iGranato et"ai1l2004l ; ICroton et al.ll2006r ). though significant 
uncertainties still affect high and low redshift measu rements 
l|Bernardi et alj|20ld ; [Marchesini et al.ll2009l , l20ld) . 

The study of number densities in other properties such 
as size, is ins tead still a t its i nfancy. The first study was 
carried out bv lCole et all ]200d). A more rec ent attempt has 
been carried out bv Ishankar et ID (|2010al ) who compared 
detailed predictions from hierarchical models with the latest 
releases of the Q(R e ) galaxy size function, finding that the 
iBower et al.l (|2006l ) model overpredicts the number of very 
compact and very large galaxies (see also lTruiillo et al]|2009l ; 
iTavlor et al]|2010l ). 

Higher redshift measurements of portions of the size 
function of spheroids are also now becoming available (e.g., 
Mancini et al.ll2009l; ISaracco et aI1l20ld : IValentinuzzi et all 



2010al ; Ivan Dokkum et alj l201Ch At higher redshifts the 



scatter in sizes at fixed stellar mass apparently seems to in- 
crease as large an d compact galaxies seem to co-exist at th e 
same epoch (e.g.. [Mancini et al.l 1201(3 ; ISaracco et alj|201ll) . 
This increase in scatter may suggest a f aster evolution in 
sizes at fixed stellar mass (|Fan et al.ll2oiol ). although trends 
between s ize and age/star formation rate, may bias this re- 



sult (e.g., i Moslch et all 1201 ll : iForster Schreiber et alJl201ll : 



ISaracco et al. 1 120111 ). 

In order to use nu mber density measure ments of size 
and velocity functions (|Cirasuolo et al.l I2005T ) at different 
redshifts, a homogeneous and well studied spectroscopic and 
photometric sample of galaxies is needed, now not yet re- 
leased. Larger statistical samples will become available in 
the near future (such as multi-wavelength optical surveys 
from the Canada French Hawaii Telescope and the Next 
Generation Virgo Cluster Survey). 

Ages, colours, metallicities, a nd other dynamical prop- 
erties (e.g.. lCappellari et al.ll2009h c an be key observables to 
probe gala xy evolut i on. For example, iBernardi et~aH (|2011al ) 
and lBernardi et all lj201ibl ). making use of the latest SDSS 
data releases, showed that all correlations wit h stellar mass 
steep e n above M Bta r > 2 x 10 n Mp) (see also iFasano et ail 
l20ld : Ivan der Wei et al] l201ll ). while relations with ve- 
locity dispersion don't. Bernardi et al. claimed, in line 
with other observational works (e.g.. lKauffmann et al.ll2004l ; 
van der Wei et all 12009) and with the results presented in 
Figure 1121 the presence of two mass scales, one at M star ~ 
3 x 1O 1O M0, below which gas dissipation controls galaxy 
formation, and a higher mass scale M star > 2 x 10 11 M B 
above which mergers dominate the evolution. The empiri- 
cal mass scales noticed by Bernardi et al. are in line with 
the characteristic masses M c and M s emphasized in Sec- 
tion [3T] Metallicity gradie nts have been recently calibrated 
by a number of groups (e .g. jFoster et al . 2009; Sp olaor et al] 
l20ld : lForbes et al.ll201 ll ) and can provide invaluable insights 
into the evolutionary patterns of early- type galaxies ( e.g., 
|Pi Matteo et al]|2009l ; II a Barbera fc de Carvalhollioog ). 



5 CONCLUSIONS 

Understanding structural evolution of spheroids has become 
one of the hottest topics in cosmology in the last years as 
it can provide invaluable insights into the true physical pro- 
cesses that regulated galaxy evolution. While angular mo- 
mentum conservation may explain many properties of discs, 
the origin of bulges is still largely debated. The situation is 
even more puzzling given that at higher redshifts galaxies 
present further disparate structural and physical properties, 
from clumpy star-forming discs, to very compact, red and 
massive galaxies. Such a complicated zoology is difficult to 
reconcile within a coherent framework of galaxy formation, 
and in fact we discussed in Section Q] that models sometimes 
propose conflicting scenarios. 

Our aim in this work was to study the predictions of 
a state-of-the-art hierarchical model of galaxy formation, 
which evolves the sizes of spheroids via mergers and disc 
instabilities, against the most recent local and high redshift 
data. To this purpose we updated the source code of the 
latest rel ease of the Munich semi-analytic galaxy formation 
model bv lGuo et al] (|201ll ). by modifying the computation 
of bulge radii exploring a variety of possibilities, and added 
the calculation of the coupled velocity dispersions. In order 
to properly compare model predictions with available data 
we made use of a large sample of early-type galaxies from 
SDSS for which bulge-to-disc decompositions have been per- 
formed both via the SDSS automated "on-the-fly" analysis 
and by applying a detailed fitting code. 

Our main results can be summarized as follows. 

• Sizes are computed in the model via energy conserva- 
tion in dissipationless mergers. Taking the model at face 
value at masses below M s < 10 M@ the model predicts 
a flattening of the size-mass relation at variance with the 
data, already at z ~ 2. 

• Following the results of hydro-simulations, we have in- 
cluded the energy dissipated in gas-rich major mergers in 
the energy budget. This modification produces progressively 
more compact remnants with decreasing stellar mass, pro- 
portionally to the fraction of cold gas in the progenitors, 
improving the match with a variety of observables. 

• We confirm and discuss evidence for two characteris- 
tic masses. One is at M c ~ 3 x 10 10 Mq , below which the 
initial bulge sizes are controlled by dissipation (higher gas 
fractions) and then evolve under mergers and disc instabil- 
ities. Galaxies above M s > 10 11 Mq, instead mainly grow 
through minor dry mergers, especially at z < 1 (Figure [T2"j) , 

• We find that the global scatter (1-ct uncertainty) in sizes 
at fixed stellar mass for galaxies with _B/T > 0.7 is com- 
parable, though systematically higher by < 40% than the 
observed one. The predicted amount of scatter for galax- 
ies with stellar mass M sta r >3x 10 10 Mq does not depend 
much on dissipation and/or amount of orbital energy in the 
merger. 

• Spheroids are predicted to be, on average, more com- 
pact at higher redshifts at fixed stellar mass. More specif- 
ically, at fixed B/T a nearly mass- independent relatively 
mild decrease in sizes is predicted, in possible disagreement 
with some observations. 

• The model predicts that environment plays a significant 
role in defining the structural properties of bulged galaxies. 
(Central) galaxies residing in denser environments are pre- 
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dieted to undergo more mergers, thus evolve faster and end 
up having larger sizes at fixed stellar mass. 

• The progenitors of the massive spheroids with M s tar ~ 
(1 - 2) x 10 11 M Q today with B/T > 0.7 are predicted to be 
compact, low starforming z ~ 2 protogalaxies with a median 
B/T ~ 0.3, with only ~ 20% remaining bulge-dominated. 
The progenitors of lower-mass spheroids with B/T > 0.7 
tend to be closer to the starforming and gas-rich proto-discs 
observed at similar redshifts. 

• Finally, we also discuss a number of related issues (in 
the text and in the Appendices), ranging from the correla- 
tions with galaxy age, with central black hole mass, with 
velocity dispersion, to the scaling relations of pseudobulges. 
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APPENDIX A: RELATION WITH AGE 

Although a large number of massive early-type galaxies 
at high redshifts have been observed to be more compact 
and with higher velocity dispersion than their local coun- 
terparts, stni_a_si£nifican^fraction shows to be already 
evolved. iMancini et al.l (|2009l ) suggested a downsizing sce- 
nario in sizes, with the most massive galaxies approaching 
the local size-mass rela tion earlier than less massive ones. 
ICappellari et all ((2009) also discussed that a large fraction 
of their galaxies at 1.4 < z < 2 with velocity dispersion 
from stacked spectrum, are co nsistent to the most dense lo- 
cal galaxies of the same mass. ISaracco et al l (|201ll l showed 
that at the average redshift of z ~ 1.5 older galaxies at fixed 
stellar mass tend to lie a factor of ~ 2 — 3 below the size- 
stellar mass relation characterizing local early-type galax- 
ies, while younger galaxies are consistent with it. The lat- 
ter type of observations introduced the concept of assembly 
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Figure Al. Predicted median size of galaxies of a given stellar mass as a function of mass-weighted age for different B/T thresholds, 
and with and without dissi pation included, as la belled. This relation for bulge-dominated galaxies is flat at fixed stellar mass, as observed 
in the local Universe (e.g,, IShankar eiT al. 20f0b). 



bias, i.e., the youngest galaxies at any epoch may be larger 
than their older coun terparts of similar stellar mass (but see 
IWhitaker et~alll2012h . 

In the local Universe the assembly bias seems instead 



to be erased, at least above M s t 



10 iu M^ and large 



BIT. IShankar fc Bernardl (|2009l ), IShankar et all (|2010b., . 
iBernardi et all |2010r ) showed from a large sample of early- 
type galaxies extracted from SDSS that old and young galax- 
ies of similar stellar ma ss share simila r sizes , i.e., the size-age 
relation is rather flat. iTruiillo et ail |201ll ). more recently, 
confirmed these results showing that the size-age relation 
is flat already at z ~ 1. When moving to lenticular galax- 
ies with lower BjT the relation gets tilt ed, with the more 
compact bulges being the oldest (e.g., Ivan der Wei et all 
120091 ; IBernardi et al.|[2O10l ). Whatever mechanism puffs up 
spheroids must be fine-tuned to allow all galaxies of a given 
stellar mass to end up on the same size-mass and size-age 
relations in the local Universe. 

Preliminary theoretical studies to interpret the size-age 
relation were perf ormed by | Shankar et al.l (|2010bl ) within 
the context of the iBower et al. ( 2006T ) hierarchical model. 
Their study showed that because older galaxies usually un- 
dergo more mergers, they naturally grow more and are able 
to "catch up" with the younger ones producing a rather flat 
size-age relation. Here we wan t to extend the preliminary 
study bv lShankar et all (|2010bh to the case of the present hi- 
erarchical model that, at variance with the Bower et al. one, 
includes dissipation and better matches the local size-mass 
relation. We will focus here only on the predicted size-age 
relation in the local Universe. At redshifts z > 1 — 1.5 the 
spread in ages becomes too small at fixed stellar mass and, 
along with its large scatter, does not provide any meaningful 
prediction to compare with the data. 

Figure IA1I shows the predicted median size of galaxies 
of a given stellar mass as a function of mass-weighted age 
for different B/T thresholds. For lower masses (left panel) 
the model predicts that in the absence of dissipation (dot- 
ted lines) more compact galaxies are also the oldest ones. 
This is mainly induced by redshift evolution in the sizes of 
the progenitors that shrinking at higher redshifts (Figure [7]) 
produce more compact early-type remnants. Including dis- 



sipation (solid lines) flattens the size-age relation out. As 
discussed above, dissipation acts in a way to reduce sizes 
proportionally to the gas fractions in the progenitors. Within 
a given mass bin, the lower mass galaxies have, on average, 
higher / gas thus suffering more dissipation than the galax- 
ies lying closer to the more massive end of the bin. Thus, in 
this model both mergers and dissipation act in a coordinated 
way to flatten out the relation. 

At higher masses (right panel) the inverse relation of 
sizes with age disappears, even in the absence of dissipation. 
In fact, including dissipation slightly steepens the relation. 
These results are all in broad agreement with what observed 
in SDSS. However, there are also some discrepancies. As an- 
ticipated above, data on lenticulars, or in general galaxies 
with lower B/T ratios, seem to foll ow a different trend, with 
older galaxies being more compact (|Bernardi et al.120101 ) . As 
seen in Figure lATl however, the model seems not to show any 
significant change in the shape of the size-age relation when 
lowering the B/T threshold (long-dashed lines; only galaxies 
with bulges mainly grown through mergers are considered 
here). This inconsistent behaviour needs to be further un- 
derstood both observationally and theoretically. 



APPENDIX B: PROPERTIES OF 
PSEUDOBULGES 

There is increasing empirical evidence that not all 
bulge s in the local Universe can be formed via mergers 
(e.g.. iKormendv fc KennicuttJ 12004 iFisher fc Drorvl 120081 ; 
iGadottil I2009F " Recent works suggest in fact that a large 
fractio n, possibly the majori ty of local bulges are pseudob- 
ulges (|Fisher fc Drorvl l201ll ). with well defined properties 
different from those of classical bulges of similar mass. Pseu- 
dobulges are usually characterized by younger stellar popu- 
lations, they are usually rotation rather than pressure sup- 
ported, have less concentrated surface brightness profiles, 
and tend to follow well distinct scaling relations in their 
global structural properties with respect to classical bulges. 

In particular, pseudobulges are found to be much more 
compact by a factor of a few at fixed stellar mass with re- 
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Figure Bl. Predicted median effective radius as a function of 
stellar mass for bulges with sizes grown mainly by mergers (long- 
dashed line with their l — o~ dispersion) , and bulges mainly grown 
via disc instability (solid lines). The classical bulges are shown 
for B/T> 0.5, while the upper and lower solid lines refer to 
pseudobulges with B/T> 0.2 and B/T> 0.5, respectively. The 
colored contours are the SDSS subsample with corresponding 
B/T> 0.5, while the triangles a nd ci rcles are pseudobulges mea- 
surements from iGadott] d2009f ) and iFisher fc Drorvl l|2010h . re- 
spectively, with no cut on B/T. 



spect to classical bulges. The open triangles and circles in 
Figure IB 1 1 are the measure d siz es of pseudobulg e s from the 
samples of iGadottil (|2009h and IFisher fc Drorvl (|2010l ). re- 
spectively. Both their stellar mass measurements have been 
converted to a c ommon initial mass fu nction using the tab- 
ulated values in iBernardi et~aD (|2010h . It is interesting to 
note that despite the completely independent and different 
methods utilized, these groups agree in finding that pseu- 
dobulges are more compact by a factor of a few. 

The model predicts a size-mass relation for pseudob- 
ulges (solid lines and filled squares) - that we here label as 
those bulges that have grown more than 50% of their final 
size via disc instability - with a slope and scatter similar to 
the one predicted for classical bulges (long-dashed line and 
open squares) but lower in normalization by a factor of a few. 
The upper and lower solid lines referring to pseudobulges 
with B/T> 0.2 and B/T> 0.5, respectivel y, are compared 
with d ata by lGado"ttil (|2009l . triangles) and IFisher fc Drorvl 
l|2010l. circles), with no cut on B/T. The predictions are in 
good agreement with the data, and even more so given that 
no fine-tuning was imposed in the model parameters. The 
model also predicts that no pseudobulges are found above 
M star > 2xl0 n M o ,a gain in line with the data. As discussed 
above, we stress that most of the pseudobulges produced in 
the model have B/T< 0.7. 

More specifically, Figur dB2l shows that the fraction of 
all pseudobulges (cyan lines) and classical (red lines) bulges 
(here all bulges are considered) are a non-trivial function 
of stellar mass. The fraction of pseudobulges tends to peak 
around Af star ~ 5x 10 10 Mq, and sharply decreases at higher 
masses, where the contribution of classical bulges domi- 
nates. Th is behaviour is in good agreement with the recent 
study by IFisher fc Drorvl (|201ll . see also Kormendy et al. 
2010), who presented an inventory of galaxy bulge types in 



Figure B2. Fraction of pseudo (cyan lines) and classical (red 
lines) bulges as a function of stellar mass for redshift 2 = (solid 
lines) and z = 1 (dotted lines). There is a tendency to have an in- 
creasing fraction of pseudobulges at lower stellar masses in broad 
agreement with the data. 



a volume-limited sample within the local 11 Mpc volume us- 
ing Spitzer and HST data (see their Figure 3) . Theoretically, 
this is expected because mergers dominate the size evolution 
above the characteristic mass of 10 11 Mq, as anticipated in 
Section 13.11 and further deve l oped in Section 14.11 On the 
other hand, IFisher fc Drorvl l|201ll ) also claim that pseu- 
dobulges tend to be the dominant class of bulges at lower 
masses, being close to 60% around M st ar ~ 10 10 M Q . The 
model instead predicts a lowering of the fraction of pseudob- 
ulges towards lower masses. More detailed observational and 
theoretical studies of pseudobulges in the low-mass regime, 
beyond the scope of the present paper, are needed to fully 
understand this discrepancy. 



APPENDIX C: THE TILT OF THE 
FUNDAMENTAL PLANE 

Early-type galaxies obey a tight scaling relation linking a, 
Mstar, and R s , the so-called fundament al plane (FP), as ex- 
pecte d from basic virial arguments fe.g.. lDiorgovski fc Davis] 
Il987l ). If the mass-to-light ratio is constant at all masses, 
i.e., L oc Mdyn, from the virial relation one would then 
simply expect L oc a 2 R e . However, the observed FP re- 
lation has small but significant departures from the theo- 
retical predictions, suggesting that Md yn /L has a non triv- 
ial dependence on mass. This "tilt" of the FP is possibly 
a consequence of stellar effects and/or progressive differ- 
ence in the inner dark matter content of t he most mas- 
sive g alaxies (e.g.. lLa Barbera et alj 120081; iTortora et alj 
20091; lLa Barbera et all |2010|; iNapolitano et al.l l20iq 



Graves fc Faberll2010l; ITortora et al.ll2012l ). 

Shankar fc Bernardil 72009 ) wed from a large sample 



of SDSS early-type galaxies that the ratio Md yri /M s tar oc 



MJ, 



with 7 ~ 0.13, i.e., it increases with increasing M s i 



in a way approximately independent of the age of the galax- 
ies. The left panel of Figure lB3l shows the predicted ratio of 
dynamical mass to stellar mass as a function stellar mass for 
galaxies of different (mass-weighted) age at z — and with 
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Figure B3. Predicted ratio of dynamical mass to stellar mass as a funct ion of Age at z = {left) and at different redshifts {right). 
The red long- dashed line is the tilt measured bv lShankar fc Bernard! l|2009h . The tilt of the fundamental plane is predicted to be rather 
independent of galaxy age at z = 0, and to flatten out at z > 0. 



B/T > 0.9. In agreement with the SDSS data, the model 
correctly predicts a positive tilt in the FP relation very sim- 
ilar to the observed one (red long-dashed line), and with 
weak dependence on age. The parametrization adopted in 
Eq. ([SJ naturally induces a tilt because lower mass galaxies 
have proportionally smaller sizes, therefore lower dark mat- 
ter fractions within Rh, thus proportionally lower velocity 
dispersions. The right panel of Figure [B3l shows that the tilt 
is predicted to decrease at higher redshifts because galax- 
ies at higher redshifts are associated to less massive dark 
matter hosts and also more compact galaxies, thus, again, 
proportionally lower dark matter fractions within Rh- The 
scatter of FP is measured to be small, thus being an impor- 
tant quantity to compare models with. However, a detailed 
study of the FP scatter is beyond the scope of this brief 
Appendix, and we will postpone it to future work. 



APPENDIX D: THE CONNECTION WITH THE 
CENTRAL BLACK HOLE: EVOLUTION IN 
SCALING RELATIONS 

Dynamical observations have revealed that super-massive 
Black Holes (BHs) are ubiquitous at the centres of most, 
if not all, local massive, bulge-dominated galaxies, with 
their mass Mbh ~ 10 6 — 10 9 M , tightly correlated with 
the mass and velocity dispersion o f the host bulge (see 
iFerrarese fc Fordll2005t IShankarlliocil . for recent reviews). It 
is therefore natural that constraining the evolution of mas- 
sive spheroidal systems contemporarily implies understand- 
ing the origin and evolution of BHs and of such tight scaling 
relations. 

The Munich SAM self-consistently follows the evolu- 
tion of BHs during the hierarchical evolution of galaxies. 
The model assumes that a fixed fraction of the cold gas 
is destabilized during merger eve nts and feeds the central 
BH (see, e.g.. iMarulli et al1l2008l for details). A parameter 
controls how much gas mass is funnelled onto the central 
BH and it is fine-tuned to reproduce the local Mbh-M ata r 
relation. The model has been found to be consistent with 
AGN luminosity functions and quasar clustering at differ- 



ent redshifts and luminosities (see details in IMarulli et al.l 
2008, 20091; iBonoli et al.ll2009l . boiol ). 



The left panel of Figure ID 1 1 shows the M^-a relation 
at different redshifts, as labelled, for a model with a com- 
puted wit h Eq. j5l). For refer ence, the grey stripe indicates 
the fit bv lTundo etail (2007) with its intrinsic scatter. The 
model produces a reasonable match to the data and we 
have checked that neglecting any halo mass dependence in a 
would have produced a flattening similar to the one observed 
for the Mstar-cr relation in Figure [9] 

Interestingly, at variance with the evolution found for 
the Mstar-o" relation, the model predicts a positive evolution 
for the Mbh-a relation, i.e., comparable or higher BH masses 
at higher redshifts at fixed velocity dispersion. This is in 
line with direct and indirect measurements of the M^-a 
at higher redshift in quasar host galaxies fe.g-. IShields et al 



2006 



Woo et al.ll2qod,l2008l;lTreu et alJIgOCffl : IShankar et al 



2009; Gasket 12009 : iBennert et al.ll201lh . Though BH host 



galaxies get more compact and thus possess higher velocity 
dispersions at higher redshifts, the redshift evolution in the 
A^bh-o" relation shows the opposite trend. 

The reason behind the opposite time behaviour of the 
two relations can be understood by looking at the right panel 
of Figure iDll which shows the predicted Mbh-M st ar relation 
for the same early-type galaxy subsample. The grey stripe 
is a linear relation of the type Mbh = 2 x 10~ 3 Mtmige with 
some scatter, i ndicative of what sugge s ted by a variety of lo - 
cal data (e.g., iMarconi fc Huntj|2003l : lHaring fc Rb3l2004 ). 
The overall agreement with the data at z = in the normal- 
ization and slope of the Mbh-Mstar relation is mostly a sim- 
ple consequence of the underlying model that assumes both 
BH accretion and star formation rate to be proportional to 
cold gas reservoir. A genuine prediction of the model is in- 
stead that the Afbh/M s tar ratio evolves at fixed stellar mass, 
possibly in a mass dependent way, i ncreasing by a factor 
of a few at higher redshifts (see also ICrotonl 120061). consis- 
tentl y with what derived by many groups fe.g.. |Pecarli et all 
2010). Thus the model predicts higher velocity dispersions 
but also more massive BHs at fixed stellar mass at higher 
redshifts, in a way to erase or even reverse the predicted evo- 
lution in the Mbh - cr relation, in good agreement with the 
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data. In other words, while most of the BH mass is accreted 
during the high-z gas-rich merger phase, the growth of the 
stellar mass and velocity dispersion of the host spheroid is 
prolonged to later times. 
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Figure Dl. Left panel: Predicted M^-a relation at different redshifts, as labelled, for a model with a dark matter mass-dependent a 
for galaxies with _B/T > 0.9. The grey stripe indicates the fit by Tundo et al. (2007) with its intrinsic scatter. Right panel: Predicted 
Mbjj-Mstar relation at the same different redshifts for galaxies with B/T > 0.9. The grey stripe shows a linear relation of the type 
Mbh = 2 X 10 — 3 jWbuige with some scatter, indicative of what suggested by a variety of local data. 
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